a = 1:1:1000;
b = 1./a;
disp(sum(b));
